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Abstract 

A new limit for the neutron electric dipole moment has been recently re- 
ported. This new limit is obtained by combining the result from a previous 
experiment with the result from a more recent experiment that has much 
worse statistical accuracy. We show that the old result has a systematic error 
possibly four times greater than the new limit, and under the circumstances, 
averaging of the old and new results is statistically invalid. The conclusion 
is that it would be more appropriate to quote two independent but mutually 
supportive limits as obtained from each experiment separately. 
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Recently, a new experimental limit for the neutron permanent electric dipole moment 
(EDM) has been reported Q with a substantial increase in sensitivity as compared to the 
previous result [[|. The increase in sensitivity results from combining the old experimental 
result with new data that has considerably worse statistical accuracy. As we discuss below, 
the averaging of the old and new results is not justified in this situation and one would 
better quote two independent results that are mutually consistent with and supportive of 
each other. It makes no sense to average the results because there are unexplained systematic 
errors in the old data that exceed the new 95% confidence limit given in [|TJ by a factor of 
at least two. 

The new quoted 90% confidence limit of \d n \ < 6.3x 10~ 26 ecm was obtained by combining 
the result of the new (1999) experiment 

d n = 1.9 ± 5.4 x 1(T 26 ecm; x 2 /i/ = 0.4 

(see Fig. 1, where v — 9) with the previous (1990) result, 

- 3 ± 5 x 10" 26 ecm (1) 

for which the error is dominated by systematic uncertainty in the data. The 95% confidence 
limit for the 1990 data is given as 

\d n \ < 12 x 10" 26 ecm. (2) 

Without attempting to correct for systematic effects, the average of the 1990 data shown 
in Fig. 1 is 

d n = -1.9 ±2.2 x 10~ 26 ecm; ^jv = 3.2 (3) 

and v = 14. Unfortunately, there is no clear way to use the discrete Rb magnetometers 
employed in the 1990 version of the experiment to correct for the magnetic systematic 
effect that is evidenced through the large x 2 jv an d bimodal character of the data; simply 
subtracting the average of the Rb magnetometer signals gives the result 

d n = -(3.4 ±2.6) x 10~ 26 ecm, j^/v = 2.2; (4) 

where the reduction in x 2 is due to the increase in the effective standard deviation due to 
the magnetometer noise, not a decrease in the scatter of the data, as is stated in 0. 
The discussion associated with Eq. (5) of M (also labelled here as Eq. (5)), 

d n = -3.4±3.9 stat ±3.1 syst x 10~ 26 ecm (5) 

does not accurately describe the source of errors in the final quoted result. In particular, 
the statistical error in Eq. (5) is obtained by multiplying the uncertainty in Eq. (4) by 
the square root of reduced % 2 , 2.6 x V2.2 = 3.9. Therefore, the statistical error given in 
Eq. (5) has a substantial systematic contribution. Because there is no physical basis for 
subtracting the average Rb magnetometer signal from the neutron EDM result (as we will 
describe), if instead we take the intrinsic neutron EDM statistical error as that given in Eq. 
(3), Eq. (5) implies a systematic error of y/3.9 2 — 2.2 2 = 3.2 (xl0~ 26 ecm); this should be 
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added in quadrature to the systematic error given in Eq. (5) yielding an overall systematic 
uncertainty of at least 



4.8 x 10~ 26 ecm (6) 

and represents a guess at how well the average systematic fluctuations averaged to zero for 
the entire 1990 data set. This systematic error, when combined with the intrinsic statistical 
error, yields the combined uncertainty given in Eq. (1). 

In fact, the final systematic uncertainty was chosen so that the 95% confidence limit 
would encompass the obvious excess scatter in the 1990 data, reflected by the large x 2 l v — 
3.2. The final 95% confidence limit was chosen to reflect this deviation in the data because it 
was accepted that there is no reason to assume the random systematic fluctuations average 
to zero; Eq. (1) above might very well be optimistic. These fluctuations are correlated with 
disassembly and reassembly of the apparatus, and therefore are possibly the result of subtle 
changes in, for example, high-voltage leakage current paths within the magnetic shields. 

We would like to take a different approach to combining the 1990 and 1999 data groups; 
this approach gives a more accurate estimate of a possible systematic error in the final 
combined number. As shown in Fig. 1, the 1990 data is better described by a bimodal dis- 
tribution. Curve A is a 1/cr 2 bimodal weighted fit to two Gaussian probability distributions 
for the 1990 data alone. Curve B shows the bimodal probability distribution for the 1990 
and 1999 combined data sets, while C is the combined data set when the standard errors of 
the 1990 data are multiplied by V3.2. It is evident that the distribution for the combined 
data set, for both B and C, is still strongly bimodal and the use of standard statistical 
techniques as in [p]] is not justified. 

This is demonstrated in Fig. 2 which gives the 68.3% and 95% confidence regions 
(A% 2 = 2.3 and A% 2 = 6.17, respectively, for two degrees of freedom @) for the two 
means determined in the fit to the bimodal distribution (a total of six parameters in the 
fit). Of particular interest is Fig. 2 A; x 2 g for the optimum fit to two Gaussian distribution 
(each Gaussian has three fit parameters so there are 25 — 6 = 19 degrees of freedom) is 8.3 
as compared to \\ 2 = 48 for a fit to a single Gaussian distribution. The reduced x 2 f°r the 
bimodal fit is Xig/19 = 0.44. However, the statistical significance is best illustrated by use 
of the F-test [Q. Adding the second Gaussian distribution is equivalent to adding three fit 
parameters, and 

xfs/19 (.44-1) v 

(the division by the new reduced x 2 1 v is valid only when it is greater than 1). The statistical 
confidence of the new bimodal description of the combined data can be calculated (see 
Appendix C-5 of ||); the probability of such a large F being due to statistical fluctuations 
is about 3%, i.e. the statistical confidence is 97%. There is no justification for artificially 
increasing the errors of the 1990 data by V3.2 because the final combined data set without 
doing so as described by the bimodal distribution has a reduced x 2 l v — 0.44. We again 
point out that even if the 1990 errors are increased by a/3T2 the full data set is still bimodal, 
as shown in Fig. 1 C and Fig. 2 B. 

The bimodal means are of opposite sign, but of remarkably close magnitudes: 
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d n i = 12.1 ± 5.0; d n2 = -11.7 ±4.5 x KT ib ecm 



(8) 



and represent the magnitude of the random systematic in the data. In Eq. (1), there is an 
implied faith that the systematic errors in the data average to zero, and therefore Eq. (1) 
must be treated with a certain caution. The Particle Data Group does indeed recommend 
scaling of the data by as done in [1]. However they justify this by the remark that this 
approach has the property that if there are two values with comparable errors separated 
by much more than their stated errors the error on the mean is increased so that it is 
approximately half the interval between the two discrepant values. Applying this criterion 
to the present case and taking the bimodal means for the 1990 data as two discrepant values 
shows a systematic error of 

d n \ d n 2 1 „ 1 n — 26 tr\\ 

= 12 x 10 ecm (9) 

which is rather larger than the weighted error when scaled by a/3.2, i.e. a factor of two 
greater than Eq. (6). This implies a 95% confidence limit of about 20 x 10~ 26 ecm, about 
a factor of two greater than Eq. (2). In addition, both x 2 plots in Fig. 2 indicate a net 
95% confidence limit of about 20 x 10~ 26 ecm. Under the circumstances, the 95% confidence 
limit given in ||, Eq. (2) above, requires additional justification. 

Furthermore, other analysis techniques support this magnitude of random systematic 
for the 1990 data. For example, when only data where the magnetometer signals were 
statistically zero was included in the average, a greater than two-sigma hint of a non-zero 
value the neutron EDM was obtained |J, 

d n = -(17 ±6) x 10~ 26 ecm. (10) 

without correcting for the magnetometers. There is no statistical or physical basis for 
subtracting the average magnetometer signal, as it has been proven by use of correlation 
techniques [|7]|| that even when the average magnetometer signal is zero, the individual 
magnetometers are still correlated with the electric field, and with the neutron EDM signal. 
This implies that the magnitude, and possibly sign, of the correlated field measured by the 
Rb magnetometers is different among the three magnetometers and the neutron storage cell. 

A better correction for a possible systematic effect was obtained by determining the 
correlation coefficient between each magnetometer and the neutron EDM signal. Subtracting 
the correlated Rb magnetometer signal made the neutron data internally more consistent; 
in this analysis, before subtracting the correlated signal, d n = —(2.8 ± 3.0) x 10~ 26 ecm, 
X 2 = 3.8, and after subtraction, d n = —(5.0 ± 3.8) x 10 -26 ecm, x 2 — 1-4. The decrease in 
X 2 in this case results from a reduction of scatter rather than an increase in the effective per 
point standard deviation due to magnetometer noise. 

Agreement between the various analyses was used to estimate the residual random sys- 
tematic in the EDM signal that did not average away, and in particular, was chosen so that 
the 95% confidence limit of the combined statistical and systematic error encompassed the 
range of systematic fluctuation, as mentioned above. Although the correlation technique is 
compelling, it is not proof that the systematic contribution has been fully accounted for. 
The implication is that the average field seen by the magnetometers does not necessarily 
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represent the field seen by the neutrons, and was the motivation for construction of the 
199 Hg co-magnetometer ||. 

To our knowledge, the correlation techniques were not applied to the new data, and 
therefore the 1999 data cannot be used to infer anything about the source of the random 
systematic in the 1990 data. In particular, we disagree with the statements associated with 
Eq. (5) of HI; the sources of the scatter in the previous data are still unknown, and ex- 
perimental evidence suggests they are associated with the magnetic fluctuations generated 
with the application of high voltage; this is the systematic of principal concern. Although 
the statements in |T| regarding magnetic fluctuations are technically correct, the "other, 
unknown, systematics" suggested in 0, if they even exist, are likely irrelevant to any dis- 
cussion. Again, the important point is whether the magnetic fluctuations that contaminated 
the data are associated with the high voltage or with other noise sources. 

The new data sheds no light on the previous result, particularly when one considers the 
vast array of changes with the new configuration of the apparatus. For example, the much 
larger ultracold neutron storage vessel in the new experiment means that a discrete magne- 
tometer placed within the shield would be relatively farther from the geometrical center of 
the neutron storage region, so the relative effect of a leakage current magnetic field could be 
quite different between the previous and present configurations of the apparatus; although 
Rb magnetometers were included in the new version of the apparatus, no presentation or 
discussion of experimental data from them is given in [|l|. Of more critical concern, in || 
the apparent 199 Hg EDM in the co-magnetometer data is not given, and is crucial toward 
understanding possible systematic effects. The graph given in Fig. 3 of [p]] might seem com- 
pelling, but tells us nothing in regard to subtle effects with application of high voltage; the 
fluctuations shown in Fig. 3 of |Q are likely primarily due to external influences. In regard 
to the old data, it was accepted that external magnetic influences on the neutron EDM 
are well-represented by the average Rb magnetometer signal, while local fields, particularly 
those associated with the high voltage supply or leakage current, are not. 

Under the circumstances, the result for the 1990 data, Eq. (5) , should be treated with a 
certain caution. The final uncertainty in this result was chosen so that the 95% confidence 
limit would encompass the magnitude of the random systematic, ~ 12 x 10 -26 ecm. Nothing 
associated with the new data would allow one to assume that the rather large random 
systematic averaged to zero in the previous data. It therefore is imprudent to average the 
old and new results. Hopefully, in the not-too-distant future, a co-magnetometer neutron 
EDM experiment will be possible with ultracold neutron density equal to or much greater 
than that obtained in connection with the 1990 version of the experiment. Only then 
will a truly new and improved limit for the neutron EDM, without concerns of systematic 
contamination of up to 12 x 10~ 26 ecm or greater, be possible. 

We thank Michael Romalis for helpful discussions. 
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FIGURES 



FIG. 1. Previous and new neutron EDM data averaged over subsequent several- week reactor 
operation cycles. The horizontal axis is arbitrary and labels the chronological order of the data 
subsets. Particulary for the 1990 data, the apparatus was dis- and re-assembled between reactor 
cycles. A: Fit to a bimodal distribution for the 1990 data. B: Fit to a bimodal distribution for 
the combined 1990 and 1999 data. C: Fit to a bimodal distribution with the 1990 data errors 
multiplied by \/3.2. 

FIG. 2. 68.3% and 95% confidence regions for the two means in the bimodal fits. A: Combined 
1990 and 1999 data. B. Combined data, 1990 data error multiplied by \/3T2. 
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